function cfdSaveUnStedyRestartDataFile(timeIrt)
%--------------------------------------------------------------------------
%  Written by LiuHaHa @ NWPU, 2022 01 
%  Contact me at: liuzhikan@mail.nwpu.edu.cn
%==========================================================================
% Routine Description:
%   This function save the unsteady resatart file
%--------------------------------------------------------------------------
global cfdFlow slvpara;
phyTemp = cfdFlow.Phy;
phyPreTemp = cfdFlow.phyPre;

filename=[slvpara.meshFile,'/Restart/','steadyRestart.dat'];
fid=fopen(filename,'w');

tempA = mod(size(phyTemp,1),3);
tempB = round((size(phyTemp,1) - tempA)/3);
for numVar=1:4
    for item=1:3:tempB*3    
        fprintf(fid,'%f %f %f \n', phyTemp(item,numVar),phyTemp(item+1,numVar),phyTemp(item+2,numVar));   
    end
    for item=tempB*3+1:size(phyTemp,1) 
        fprintf(fid,'%f ', phyTemp(item,numVar));  
    end
    fprintf(fid,'\n');  
end

for numVar=1:4
    for item=1:3:tempB*3    
        fprintf(fid,'%f %f %f \n', phyPreTemp(item,numVar),phyPreTemp(item+1,numVar),phyPreTemp(item+2,numVar));   
    end
    for item=tempB*3+1:size(phyPreTemp,1) 
        fprintf(fid,'%f ', phyPreTemp(item,numVar));  
    end
    fprintf(fid,'\n');  
end   

fprintf(fid,'%f %f\n', cfdFlow.vel(1),cfdFlow.vel(2)); 
fprintf(fid,'%f \n', cfdFlow.PhyTotalTime); 
fprintf(fid,'%d \n', timeIrt); 
fprintf(fid,'%d \n', slvpara.flx_order); 

fclose(fid);      

end
       